A Measurement of the Branching Ratio and Asymmetry of the Decay E° — > S°7 
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We have studied the rare weak radiative hyperon decay 
H° — > E°7 in the KTeV experiment at Fermilab. We have 
identified 4045 signal events over a background of 804 events. 
The dominant H° — » A7r° decay, which was used for normal- 
ization, is the only important background source. An analysis 
of the acceptance of both modes yields a branching ratio of 
Bi?(H° -> E a j)/BR(~° -> Att°) = (3.34±0.05±0.09) x 10 -3 . 
By analyzing the final state decay distributions, we have also 
determined that the E° emission asymmetry parameter for 
this decay is ass = —0.63 ± 0.09. 

PACS numbers: 13.30.Eg, 14.20.Jn 

We report here on the Fermilab KTeV experiment's anal- 
ysis of the decay S° — > £"7. This is an example of one 
of the most intriguing classes of baryon decays- weak 
radiative hyperon decays (WRHD) . These decays are ex- 
perimentally quite accessible, having branching ratios on 
the order of 10 -3 0. The two most significant variables 
that can be measured in these decays are the branching 
ratio (BR), and the asymmetry of the baryon emission 



with respect to the initial spin axis (a). However, despite 
this experimental simplicity, there is no known theoreti- 
cal framework for explaining these decays. Some of the 
difficulty in explaining WRHD is because predictions us- 
ing the quark model do not match predictions from an 
analysis at the hadron level 

Hara proved in 1964 that a = in the SU(3) limit for the 
S + and S~ WRHD, assuming only CP invariance and U 
spin symmetry (s and d quark exchange symmetry) ||. 
An estimate, based on single-quark s — > d transitions and 
that takes into account SU(3) breaking predicts a modest 
positive asymmetry, with the magnitude depending on 
the assignment of quark masses [|J. However, the only 
asymmetry of any WRHD that has been measured accu- 
rately, the one for S + — * pj, has been found to be highly 
negative, —0.76 ± .08 ||. The corresponding branching 
ratio cannot be understood as due to single-quark tran- 
sition. A recent paper suggests that this large negative 
asymmetry could result from intermediate resonances in 
the decay, and gives a prediction of a positive asymmetry 
for the S° -> S°7 mode [§. 

To constrain the theoretical models, it is vital to measure 
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accurately the parameters of WRHD for hyperons other 
than the E + . Reference (7) states that whether Hara's 
theorem can be incorporated into a succesful theoretical 
model depends on the experimental results from the two 
S° weak radiative decays. 

The decay S° — ► £°7 has been previously observed in 
two separate experiments [^[^). Both experiments have 
sample sizes of less than 100 events. The most accu- 
rate branching ratio is reported as (3.6 ±0.4) x 10~ 3 ||. 
This experiment also reported an asymmetry measure- 
ment, but it is incorrect because the depolarization in 
the decay S° — > A7 was not properly taken into account 
as discussed below. 

Many of the details of the KTeV beamline, experiment 
and hyperon trigger can be found in a previous publi- 
cation |l(J. We reemphasize here some of the details 
particularly relevant to this measurement. The beam- 
line for KTeV was designed to deliver two square high 
intensity beams of Kl particles for CP violation studies. 
The detector is situated 94 meters from this target. The 
sweeping magnets in the beamline were designed and op- 
erated so that the integrated magnetic field delivered S° 
hyperons polarized (with about 10% polarization) in the 
positive or negative vertical direction. We reversed one 
of the magnets regularly so that the net polarization was 
zero for the data discussed here. 

The detector consists of a 64 m long vacuum decay vessel 
followed by a spectrometer. The decay vessel contains 
several lead/scintillator vetos with square apertures in 
their center. These veto events where decay particles 
leave the sensitive area of the detector. The vacuum ves- 
sel is followed by 4 drift chambers, two situated on either 
side of an analysis magnet. This magnetic spectrome- 
ter gives a momentum resolution for charged particles of 
cr(P)/P = 0.38% © P x 0.016%, where P is in GeV/c. 
Photon veto detectors surround each drift chamber sta- 
tion. 

Downstream of the last drift chamber is a set of 9 transi- 
tion radiation detectors (TRD). A plane of vertically ori- 
ented trigger hodoscopes is situated after the last TRD 
chamber. This plane is followed by the electromagnetic 
calorimeter which is composed of an array of 3100 crys- 
tals of pure Csl, 27 radiation lengths deep. The energy 
resolution of the calorimeter for electrons and photons 
is approximately a(E)/E w 0.45% © 2%/y r E and the 
position resolution is approximately 1 mm. 
There are two beam holes in the TRD radiators, the trig- 
ger counters and the calorimeter. After the calorimeter is 
a 10 cm thick lead wall and then a hodoscope array (HA), 
with a single beam hole encompassing both beams, which 
provides a veto for hadronic showers. This is followed by 
three steel walls and two hodoscope arrays, which pro- 
vide a veto for muons. Two small scintillation counters 
are situated in the beam hole of the first steel wall to pro- 
vide a trigger for highly forward going charged particles. 
The final decay products of the S° — > £°7 signal are a 



proton in one of the beam holes and a ir~ and two pho- 
tons hitting the calorimeter. These are the same final 
decay products as for the normalization mode S° — > Air . 
The trigger for the experiment is built on a standard 
multi-level architecture, with the Level 1 incorporating 
fast triggering elements, Level 2 being a somewhat slower 
and more detailed hardware trigger and Level 3 being 
a real time processing filter to select for specific decay 
modes. The main hyperon trigger was designed to mea- 
sure the previously unobserved beta decay of the H°. 
Further details can be found in the paper reporting its 
discovery p0|. Since this trigger vetoed on a signal of 
more than 2.5 MIP's in the HA array, it rejected 58% of 
the S° — > E°7 signal and S° — > Air normalization mode 
because the final state ir~ often showers in the Csl or 
lead. This introduces no bias in the measurement of the 
branching ratio, however, because of the identical final 
state between signal and normalization mode. 
We established a set of offline selection cuts to accept 
both modes and another set of cuts to distinguish be- 
tween the two. The cuts in common include the following. 
There are 2 charged tracks, with the higher momentum 
track positively charged and with momentum between 85 
and 600 GeV/c. The lower momentum track is negative, 
with a momentum between 5 and 150 GeV/c and has an 
energy in the Csl calorimeter less than 90% of its mo- 
mentum. The two tracks must have a momentum ratio 
greater than 3.5 and, assuming a proton and pion iden- 
tity, must combine to form an invariant mass within 15 
MeV/c 2 of the A mass. Each of the two calorimeter clus- 
ters is required to have at least 2 GeV of energy and no 
charged track pointing at it. The sum of the two energies 
must be at least 16 GeV. We also require that the two 
photons do not reconstruct as a tt° within 3m along the 
beamline direction of the charged track vertex. 
We reconstruct the A flight path from the reconstructed 
momentum of the proton and pion tracks. To find the 
decay vertex of the parent H°, we numerically vary the 
distance along the beamline, or z. For a S° vertex at a 
given z, we calculate the transverse coordinates by ex- 
trapolating the A flight path back to that location. We 
choose the value of z that minimizes the difference be- 
tween the reconstructed invariant mass of all the final 
state particles and the S° mass. We discard events with 
a decay vertex outside of the vacuum decay region, which 
we define as Z=95-155 m. Events are also discarded if the 
daughter A vertex falls upstream of the parent S° vertex. 
We also require that the total reconstructed momentum 
(A plus two photons) transverse to the flight path of the 
S° parent, or pj,, be less than .0005 GeV 2 /c 2 . In addi- 
tion, we reconstruct the event as a decay of a A from the 
primary target and as a S° — * A7 decay and discard the 
event if it is consistent with either. 

The distinguishing characteristics between the signal and 
normalization decay modes reside in the invariant masses 
between pairs of the final state particles. The signal 
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mode will have one of the A7 combinations reconstruct 
as a E°, while the normalization mode will have a 77 re- 
construction near the n° mass. Figure 1 shows these 
invariant mass combinations for the data. A cut of 
|m 77 — m^o I < 20 Mev/c 2 is made in this figure to make 
the signal mode visible. This cut gives the optimum value 
for signal to background. The lower energy photon is cho- 
sen to be paired with the A. The presence of a S° mass 
band can be seen in this plot along with the band arising 
from the incorrect pairing of the photon. 
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FIG. 1 . A plot of A — 7 mass versus 7 — 7 mass for all data. 
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FIG. 2. A plot of the A — 7 invariant mass for all data 
events passing selection cuts. Superimposed is normal mode 
Monte Carlo events that pass the selection cuts. This latter 
distribution is scaled to fit the data. 



Figure 2 shows the final invariant A7 mass plot (a pro- 
jection of Figure 1), with the photon chosen to give the 
lowest mass. The background under the S° mass peak 
can be explained entirely by the normalization mode de- 
cays S° — > A/T° that fall outside of our 77 mass cut. Wc 
fit the distribution in Figure 2 above a mass value of 1.16 
GeV/c 2 to a combination of Monte Carlo simulations of 
the S° — > S°7 signal and S° — » At: background. The 
relative level of the background has been allowed to float 
in the fit since the simulation is not guaranteed to be ac- 
curate in describing the population at the extreme tails 
of the 7r° mass peak. The fit level is about 25% higher 
than indicated by the number of well reconstructed nor- 
malization mode decays. 

As can be seen in Figure 2, the simulated background 
describes the data very well. We subtract the simu- 
lated background distribution from the data and count 
the number of events in the £° mass region of 1185- 
1201 MeV/c 2 . For our data, this is 4045 events, which 
gives a statistical uncertainty of ±64 events. The num- 
ber of background events in this region as estimated by 
the background fit is 804 events, giving a signal to back- 
ground ratio of 5:1. The number of normalization events 
in our data is 1377642. There is a negligible amount 
of background for the normalization mode. Using the 
simulation of the detector, we determine that the ratio 
of acceptances for the signal decay mode to the nor- 
malization decay mode is 0.88. The acceptance does 
not depend significantly on the asymmetry parameter 
used. This results in a branching ratio measurement of 
BR(E° -> S°7)/Bi?(S° -> Att ) = (3.34 ± 0.05) x 10~ 3 
where the error quoted is purely statistical. 
The two largest systematic effects on the calculation of 
the branching ratio are background and acceptance un- 
certainties. The statistical error on the amount of back- 
ground is 28. We varied the region around the signal 
mass peak in which to fit the simulated background and 
obtained an uncertainty of 40 events in the signal re- 
gion due to the level of the fit. Additionally, looking 
at comparisons between event variables for those back- 
ground events in the E° mass peak sidebands, we can 
see small discrepencies between the simulation and the 
data. We conservatively estimate that at most 60 back- 
ground events could be unaccounted for under the signal 
peak. In all, we assign an uncertainty of 80 events to 
the amount of background, giving an uncertainty in the 
branching ratio of ±0.07 x 10~ 3 . 

The uncertainty in the acceptances is primarily due to 
the fact that we do not simulate hadronic showering in 
the HA detector and thus do not account for this effect 
at the trigger level. The effect of the hadronic veto on 
the absolute level of our signal is substantial. However, 
because the normalization mode has the identical final 
state, this effect cancels out in the branching ratio calcu- 
lation. We find that the kinematic distributions relevant 
to the HA detector (such as pion and proton position 
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and momentum) are very similar between the signal and 
normalization decay modes. We tested a simple model 
of hadronic vetoing in the Monte Carlo containing both 
position and momentum dependence for the tt~ . We ob- 
served that the ratio of acceptances between signal and 
normalization did not change at all. We assign a value 
of ±0.05 x 10 -3 branching ratio uncertainty due to the 
hadronic veto. 

We explored the sensitivity to analysis cuts by varying 
the three most relevant variables: p\, the z range for the 
S° and the separation of the S° and A vertices. The 
variation seen from these cuts matches the level of un- 
certainty from the other sources quoted above. We also 
broke the data into various subsets (left vs right beam, 
and run number) and find the individual results statisti- 
cally consistent. 

The final systematic uncertainty in the branch- 
ing ratio from the above sources is ±0.09 x 10~ 3 . 
This then gives the final branching ratio result of 
BR{E° -► S° 7 )/ J Bi?(S -► Att ) = (3.34±0.05±0.09) x 

10- 3 . 

A value of the asymmetry parameter, owe , has also been 
obtained from our data. This parameter is determined 
most directly from the up-down asymmetry in the de- 
cay distribution of polarized H°'s. However, our low de- 
gree of polarization (w 10%) limits the power of this ap- 
proach. For an unpolarized initial S° state, as in our to- 
tal data sample, angular momentum considerations dic- 
tate that the longitudinal polarization of the decay E° 
is just — ose In the subsequent electromagnetic de- 
cay, the A retains the component of the parent S° po- 
larization along the A direction of motion, but with op- 
posite sign. The final A — > p7r~ weak decay asymmetry 
(«Ap = 0.642 pi) then analyzes the polarization. The re- 
sulting proton angular distribution is thus proportional 
to [1 + as£aApCos(#£A)cos(0Ap)], where the angles are 
determined in the rest frame of each decaying particle. 
If one ignores the decay S° — > A7 and simply averages 
over all directions of the A emission, then the net asym- 
metry seen will be diluted by a factor of -1/3, assuming 
that experimental acceptance does not affect the decay 
asymmetrically. It should be noted that the previous 
measurement of asymmetry in the decay S° — > S°7 did 
not take into account this dilution [pi. 
In our analysis, we consider the 2-dimensional distribu- 
tion of the cosines described above. We measure the rest 
frame angles for each event and then perform a \ 2 com- 
parison between the distribution seen in data with that 
for Monte Carlo, using different values of the asymmetry 
parameter in the generation. We used a 5 bin by 5 bin 
construction. A row of 5 bins is virtually empty due to 
the cut made on the tt° mass so we do not use these bins 
in the x 2 calculation. We made 11 separate Monte Carlo 
samples using different values of the asymmetry param- 
eter, ranging from owe = 0.0 to asx: = —1.0. The x 2 
between the background subtracted data and the Monte 



Carlo distributions was calculated for each of these cases 
and the result is shown in Figure 3. A fit to this curve 
with a parabola gives a minimum at ass = —0.63, with 
a corresponding value of x 2 of 27.7 for 18 degrees of free- 
dom. The 1 unit variation of x 2 gives an error on this 
measurement of ±0.08. 

Several tests were done to see if there are systematic ef- 
fects on this result for the asymmetry parameter. To 
test the effect of background on the result for asymmetry, 
we repeated the analysis described above without back- 
ground subtraction. The result is ol&e. = —0.66 ± 0.08, 
which does not vary from the background subtracted 
value. The main difference is that the % 2 /DOF at the 
minimum increases from 1.5 to 2.8. To test whether 
the beam hole region had an effect on the calculation 
of the asymmetry value, we applied the additional re- 
quirement that the pion does not hit the region of the 
calorimeter between the beam holes. We then repeated 
the x 2 analysis. The result is ckh£ = —0.67 ± 0.08 with 
the x 2 /DOF at the minimum of 1.5. We also calculated 
the x 2 fit for comparisons with the Monte Carlo data 
without HA simulation and found no significant differ- 
ence. We estimate that the systematic error on the asym- 
metry measurement is ±0.05, yielding our final result, 
«hs = -0.63 ± 0.08 ± 0.05. We independently verified 
the sign and magnitude of this asymmetry by directly 
comparing our two oppositely polarized samples of data. 
This result is only the second measurement making an ac- 
curate determination of a WRHD asymmetry. The value 
of the decay asymmetry is highly negative, similar to that 
seen in the decay E + — * pj. 



a = -0.63 
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FIG. 3. The plot of \ /DOF comparison of background 
subtracted data cosine distributions to the Monte Carlo. The 
2nd degree polynomial fit is shown, with the minimum at 
a = —0.63 identified. 
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